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Abstract: A hypergeometric type equation satisfying certain conditions defines either a 
finite or an infinite system of orthogonal polynomials. The associated special functions are 
eigenfunctions of some shape invariant operators. These operators can be analysed together and 
the mathematical formalism we use can be extended in order to define other shape invariant 
operators. All the considered shape invariant operators are directly related to Schrodinger type 
equations. 
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1 Introduction 

Many problems in quantum mechanics and mathematical physics lead to equations of the 
type 



where a{s) and r(s) are polynomials of at most second and first degree, respectively, and 
A is a constant. These equations are usually called equations of hypergeometric type [7j, 
and each of them can be reduced to the self-adjoint form 
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by choosing a function g such that (erg)' = tq. 
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The equation (jT]) is usually considered on an interval (a, b), chosen such that 

cr(s) > for all s G (a, b) 

g{s) > for all s G (a, b) (3) 
lims^acr{s)g{s) = lims^b<y{s)g{s) = 0. 

Since the form of the equation (0) is invariant under a change of variable s cs + c?, it is 
sufficient to analyse the cases presented in table 1. Some restrictions are to be imposed 
to a, (3 in order the interval (a, b) to exist. The equation (jT} defines an infinite sequence 
of orthogonal polynomials in the case a{s) G {1, s, 1 — s^}, and a finite one in the case 
a{s) G {s^-l, s\ s^ + l}. 



a{s) 


t{s) 


q{s) 


a, (3 


(a, 6) 


1 


as + (3 




a < 


M 


s 


as + j3 




a < 0, /3 > 


(0,oo) 




as + j3 


(1 + 5)-(a-/3)/2-l(i _ ^^^)-(a+/3)/2-l 


a < /3 < —a 


(-1,1) 




as + j3 


(S + l)(a-/3)/2-l(^ _ ;l)(q+/3)/2-1 


-/? < a < 


(l,oo) 




as + (3 


^a-2g-/3/s 


a < 0, /? > 


(0,oo) 


+ 1 


as + (3 


_j_ ^2^o/2— arctan s 


a < 


R 



Table 1 The main cases 



The literature discussing special function theory and its application to mathematical 
and theoretical physics is vast, and there are a multitude of different conventions con- 
cerning the definition of functions. A unified approach is not possible without a unified 
definition for the associated special functions. In this paper we define them as 

^iAs)={^Ms)y^Us) (4) 

where are the orthogonal polynomials defined by equation ((H). The table 1 allows 
one to pass in each case from our parameters a, /? to the parameters used in different 
approach. 

In 121 El we presented a systematic study of the Schrodinger equations exactly solvable 
in terms of associated special functions. In the present paper, based on the factorization 
method PI El and certain results of Jafarizadeh and Fakhri P , we extend our unified 
formalism by adding other shape invariant operators. 
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2 Orthogonal polynomials 

Let r(s) = as + /3 be a fixed polynomial, and let 



a" is),,, a" 



2 l{l-l)~r\s)l=--l{l-l)-al (5) 

for any / G N. It is well-known ^ that for A = Aj, the equation admits a polynomial 
solution $i = of at most / degree 

a{s)<^'l + r{s)^ + \i^i = Q. (6) 

If the degree of the polynomial $/ is / then it satisfies the Rodrigues formula [Zj 

where Bi is a constant. Based on the relation 

{ (5 G M I lim,^aa(s)^(s)/ = lim.^b (t(s)^(s)s^ = } 

[0,oo) if (t(s) G {1, s, l-s2} (8) 
[0,-a) if G {s^ - 1, s^, s^ + l} 

one can prove [HIIZI that the system of polynomials {$; | I < A}, where 

oo for cr(s) G {1, s, 1 — s^} 

A=<; (9) 

^ for (t(s) G s2 + 1} 

is orthogonal with weight function q{s) in (a, h). This means that equation defines an 
infinite sequence of orthogonal polynomials 

^"0, $1, $2, ••• 

in the case a{s) G {1, s, 1 — s^}, and a finite one 
$0, '^'i, 

with L = max{/ G N | / < (1 - q;)/2} in the case a{s) G {s^ - 1, s^, + 1}. 
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The polynomials ^^"'^^ can be expressed (up to a multiplicative constant) in terms of 
the classical orthogonal polynomials as 



y 2 



-2a 



,(-(a+/3)/2-l, {-a+/3)/2-l). 



,({a-/3)/2-l, (a+/3)/2-l), 



I 



J l-a-2l (13 



-s 



(!) 



in the case a{s 

in the case cr(s 

in the case cr(s 

in the case cr(s 

in the case a(s 



.,p((„+i;3)/2-l, (a-i;3)/2-l)^.^^ ^j^^ ^^^^ 



1 

S 



(10) 



where H;, Lf and pj^''^-* are the Hermite, Laguerre and Jacobi polynomials, respectively. 
The relation (fTIHl does not have a very simple form. In certain cases we have to consider 
the classical polynomials outside the interval where they are orthogonal or for complex 
values of parameters. 



3 Associated special functions. Shape invariant operators 

Let / G N, / < A, and let m G {0, 1, /}. The functions 

d'^ 



Ms) 



where 



are called the associated special functions. If we differentiate | 
multiply the obtained relation by k"^{s) then we get the equation 

where is the differential operator 



V(s) (11) 
3) m times and then 



Ms) 



ris)i + 



d_ I m(m-2) ((T'{s)f 
(is ~^ 4 fr(s) 



I mT(s) cr'(s) I / n\ II I \ I { \ 

+ ~ - 2)0- (s) - mr'(s). 



(12) 



(13) 



For each m < A, the special functions with m < I < A are orthogonal with 
respect to the scalar product 



if, 9) 



f{s) g{s)g{s)ds 



(14) 



and the functions corresponding to consecutive values of m are related through the rais- 
ing/lowering operators |21 E] 



Am = «(s)^ - mK'{s) 

At^ = -n{s)i-^^^-im-l)n'is) 



(15) 
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namely, 

{0 for / = m 

^i,m+i for m < / < A (16) 
^+$«,m+i = (A«-A^)<l>«,^ for 0<m</<A. 
Up to a multiplicative constant 

kHs) for m = I 

^Us) = { y ^ (17) 

An -^^=^kHs) for m<l 



A;— Am A;— Am + 1 "' A;— A;_i 

and the operators are shape invariant [21 E] 

Hm+1 ~ Am = A^A^ H^A^ = A^i^m+l- 

The functions 



(18) 



^/,m = $«,m/||$/,m|| (19) 



where ||/|| = a/ (/, /) are the normalized associated special functions. Since |21IS] 

||$i,m+l|| = VAz-Am||$Z,m|| (20) 

they satisfy the relations 

for Z = m 



m 'rLm 



VA; - Am 0z,m+i for m < / < A 



(21) 



0i,m+l = VA/ - Am 0/,m for < 771 < / < A 

= — /-+^ ... /''-^ (/>n. 

4 Application to Schrodinger type equations 

It is well-known j5| that the equations ifm$z,m = Ai$;^m are directly related to certain 
Schrodinger type equations. If in equation satisfied by $i^m 



(22) 



-a(.)^$,m(.) - r{s)£<l>Us) + 

+ - ^^(^ - 2)a"(s) - mr'(s)]$,,m(s) = A,$,,m(s) 

we pass to a new variable x = x{s) and a new function ^^^^(x) such that 

dx 

^ = as) <^l,m{s) = vis) ^l,m{<s)) (23) 
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^(s) 7^ and 77(5) 7^ for any s G (a, b), then we get the equation 
where 

T/ re^l — "^("^-^) ( 

a{s) ^ 2 tT(s) 



(24) 



(25) 



and the dot sign means derivative with respect to x. For C,{s) and r]{s) satisfying the 
conditions 

ais) e{s) = 1 ais)eis) + 2a(.) ^(5) ^ + r(.)e(.) = 0. (26) 

r]{s) 

which lead to 

m = ±^ vis) = ^--^ (27) 

(up to a multiphcative constant), the equation becomes 

- ^z,,,(x(s)) + VUs)'^i,m{x{s)) = Xi ^«,™(x(s)). (28) 
Denoting by s{x) the inverse of the function (a, b) — > (a', 6') : s 1-^ x{s) we get 

— = ±k{s{x)) *«,m(a;) = Vk(s(x))^(s(x))$z,™(s(x)). (29) 

The equation is satisfied for any s G (a, 6) if and only if 

- ^/,m(a;) + Kn(s(x))^',,„(x) = A, ^i,^(x) for any x G (a, 6') (30) 
that is, if and only if \E'i ^(x) is an eigenfunction of the Schrodinger type operator 

nm = —^ + Vm{x) where V^x) = V{s{x)). (31) 
For each m < A the functions \E'/,m(a;) with m < I < A are orthogonal 

X!','^«,m(a;)^fc,m(x)dx = /^^ $i,„(s(x))$fc,^(s(a;))^(s(x)) |£| 
for k I, and satisfy the relations 



^m*«,m(a;) 



for / = m 

for m < / < A (32) 
A+'^i,m+iix) = {\i - A^)^i,^(x) 



Nicolae Cotfas / Central European Science Journals () 1-13 



7 



where 

are the operators corresponding to Am and A^. Particulary, we have [S] 

n m+l AmAjyi^ 'HmA,^ A^'hCm+1- 

and 



for each m G {0, 1, Z — 1}. 



(33) 



(34) 



v^,.(x) = vi;,,(,:) (35) 



Theorem 4.1. If the change of variable s = s{x) is such that ds/dx = ±/t(s(x)) then 



dx 



Ai = T-^ + Wmix) 



and 



VUx)=WMTWm{x) + Xr 

where Wm{x) is the superpotential P 



^m,m ('^) 



2/t(s(x)) V 2y dS ^m,m{x) 

Proof. From (erg)' = rg and ds/dx = ±k(s(x)) we get 

g' T ^k' d ^ 1 d 

g K,'^ K, ds /t(s(x)) dx 

whence (jHE))- Since Am^m,m = 0, from (jSH), and (jHE)) we obtain 

i^m,m ~l~ ^'^m(''')^m,m ^m,m ~l~ (^^m('^) -^r7i)^m,m 0- 

The functions 



tpl,m{x) = a/ k{s{x)) gis{x)) (pl,mis{x)). 

corresponding to 0; m are normahzed 

-6' 



/ \i^k,m{x)\ 
J a' 



'dx 



\4>k,m{s{x))\ g{s{x)) 



ds 


dx = 


dx 


J a 



(36) 



(37) 



(3^ 



(39) 



'>k,m{s)\'^ g{s)ds = 1 
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and satisfy the relations 

for I = m 



Lm 



- Am for m < / < A 

ipi,m+i = \/k - Am i^i,m for < m < / < A ^'^^^ 

Particular cases [Ulllini- Let am = —(2m + a — l)/2, = (2m — a — l)/2. 

(1) Shifted oscillator 

In the case a{s) = 1, the change of variable M — > M : x i— »• s{x) = x leads to 



Wmix] 



ax+P 



(41) 



Vm(x) = ^^ + f + Am 

where Am = —am. 

(2) Three-dimensional oscillator 

In the case cr(s) = s, the change of variable (0, oo) — > (0, oo) : x ^ s{x) = 
leads to 

l^m(x) = -fx-(/5 + m-i)i ^^^^ 
Vm(x) = + (/5 + m - i) (/? + m - I) + f (/? + m) + Am 

where Am = —am. 

(3) Pdschl-Teller type potential 

In the case (j(s) = 1 — s^, the change of variable (0, tt) — > (—1, 1) : x ^— s(x) = cosx 
leads to 

Wmix) = a' cotanx - | cosec x = cotan f - ^%^tanf 

Vm(a;) = ^ttm^ ~ '^m + "t) cosec — (2a^ — 1) | cotanx cosecx — a^^ + Am 

where Am = m(m — a — 1). 

(4) Generalized Pdschl-Teller potential 

In the case (t(s) = s^— 1, the change of variable (0, oo) — > (1, oo) : x i-^ s(x) = coshx 
leads to 

Wmix) = am cotanh x — | cosechx 

Vm(a;) = ( am + t^m + ^ ) cosech ^x— (2am + l)f cotanh x cosech x+a^ + Am 
where Am = —m{m + a — 1). 
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(5) Morse type potential 

In the case a{s) = s^, the change of variable 



(0, oo) : X I— > s{x) = leads to 

(45) 



V^(x) = ^e-2- - (2a^ + l)fe-- + < + A, 



where = —m{m + a — 1). 
(6) Scarf hyperbolic type potential 

In the case (j{s) = s^+l, the change of variable 



Wm{x) = a^tanhx — | sechx 



: X ^ s{x) = sinhx leads to 



(46) 



Vm(a;) = (^—C(m ~ (^m + x) ^cch — {2am + l)f tanhx sechx + + A^- 
where Am = —m{m + a — 1). 

5 Other shape invariant operators 

In this section we restrict us jSj to the particular non-trivial cases when a and (3 are such 
that there exists A; G M with q{s) = cr'^(s) (see table 2). 







g{s) 


k 




(a, 6) 


s 


p 




a- 


1 


(0,oo) 




as 


(1 _ 52)-a/2-l 


a 
2 


- 1 


(-1,1) 




as 


(s2 _ i)a/2-l 


a 
2 


1 


(l,oo) 




as 


^a/2-1 


a 
2 


1 


(0,oo) 




as 


(s2 + i)a/2-l 


a 
2 


1 


(— oo, cxd) 



Table 2 The cases when g{s) = a''{s) 
From {ag)' = tq we get r(s) = (A; + l)cr'(s) = 2{k + 1)k(s) k'(s), and 

Am = K{s)i - mK'is) A+ = -k{s)£ -i2k + m + 1)k'{s) 
Hm = -t^^s)-^ - 2{k + 1)k{s) K'{s)j- - m{m + 2k)K{s) k"{s) 



(47) 



= -m{2k + m+ 1)^ = -m{2k + m + l)[/t'^(s) + /t(s) /t"(s))]. 
Theorem 5.1. If a and j3 are such that f)(s) = a^(s) then for any 7 G M the operators 



4—4-1- 



7 



2m + 2A; + 1 



m m ' 



7 



2m + 2A; + 1 



(48) 
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satisfy for m < A — 1 with 2m + 2/c + 1 7^ the relations 

A A+ — H . — \ n A+ — A+ H 

where 



(49) 



Proof. Since A^A^ = — Xm and A^A'^ = Hm+i ~ "we obtain 

(A+ + e){A^ + £) = Hm- X„,- e{2m + 2k + 1)k'(s) + 
(Am + e){Al + e) = Hm+i - \m - e(2m + 2A; + l)/s:'(s) + 

for any constant £. If we choose e = l/(2m + 2k + 1) then we get 



HmA^ (^A^Ajji -\- XjYi)A^ A^(^Afj-iA^ -\- Am) A^^^H^^i 

AjnHfji Auii^A^Ayn ~\~ Am) ^AmA^ -\- Am)Am -^m+l^m' 

Theorem 5.2. If 0<m</<A and if ^i^i satisfies the relation Ai^ii = then 

A+ 4+ 4+ 4 + 

$^ „^ = ^An_ , ... , ^^-^ . (51) 

A; — Am A; — Am+1 A/ — A/_2 A; — A/_i 

is an eigenfunction of Hm 

Hra^l^ra = Xl^l,m (52) 

and 

{0 if m = / 

<|),,m+i if m < / (53) 

i+l'/.m+l = (A; - Am)$/,m. 

Proof. The definition ()5H) of $/ m can be re-written as 

hm = ~ ^\ (54) 
Ai — Am 



and Hi(^i^i = {AlAi + A;)$i,« = A/^;,;. The relation Hm^i^rn = Xi^i,m follows by recurrence 

n A+ A+ H , 

Xl — Am Ai — Am 
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From the relation we get 

A; — Am Ai — Am 

If in equation (j52|) we pass to a new variable x = x{s) such that dx/ds = ±1/k{s) 
and to the new functions 

^i,mi^) = ^/K{s{x)) g{s{x))^i,mis{x)). (55) 

then we get the Schrodinger type equation 

nm^i,m = ~Xl^l,^ (56) 

where 

~ d'^ ~ ~ dK 

= + Vm(x) and Vm{x) = Vm{x) - 7-^(-s(a;)). (57) 



The operators 
corresponding to A^ and satisfy the relations 



Am'^l^mix) = 

A+-^l,rr,+l{x) = {\i - \rn)'^l,m{x) 



(5^ 



if m = / 

^;,m+i if m < / (59) 



(60) 



and 

"^m Am Ar^Am Am'Hm 'Hm-\-lAm 

T~^m+1 Am AmA^ HmAjy^ A^'M.m+l' 

If the change of variable s = s{x) is such that ds/dx = ±k(s(x)) then 

A^ = ±^ + WUx) A+=T^ + WUx) (61) 

V™(x) = W^X^) T ^U^) + Am = + ^™ (62) 

Vm,m(a;) 

where iym(a;) is the superpotential P 

yVm\x) = , , — m -;-{s(x)) H ; = =F-3 — ■ . 63 

Ms{x)) V 2^^^^^^^ 2m + 2k + l ^vl/m.m(x) ^ ' 



12 



Nicolae Cotfas / Central European Science Journals () 1-13 



Particular cases jHIHEl- Let = —{2m + a — l)/2, a'^ = (2m — a — l)/2. 

(1) Coulomb type potential 

In the case cr(s) = s, the change of variable (0, oo) — > (0, oo) : x i-^ s{x) = x'^/i 
leads to 

V^(x) = (/3 + m-l) (/3 + m-|)^-7i (64) 



(2m+2/3-l)2 • 

(2) Trigonometric Rosen-Morse type potential 

In the case a{s) = 1 — s^, the change of variable (0, vr) — > (—1, 1) : x ^— s{x) = cosx 
leads to 



Wmix) = a' cotan x + 



7 



2m— a— 1 



= {cx'm^ — ct'm) coscc + 7 cotan X — a'^^ + m(m — a — 1) (65) 
= m{m - a - 1) - 

(3) Eckart type potential 

In the case a{s) = s^— 1, the change of variable (0, oo) — > (1, oo) : x t-^ s{x) = coshx 
leads to 

Wm{x) = am COtanh X + ^m+a-l 

Vm(x) = (a^ + am) cosech — 7 cotanhx + — m{m — a — 1) (66) 
= -m{m - a - 1) - (^^^^^^p ■ 

(4) Hyperbolic Rosen-Morse type potential 

In the case (j(s) = s^ + l, the change of variable M — : x i— > s(x) = sinhx leads to 

W^m(a;) = Omtanhx + 



^mix) = — (a^ + am) sech^x — 7tanhx + — m{m — a — 1) (67) 
~Xm = -m{m - a - 1) - (2,jli)2 - 

6 Concluding remarks 

Most of the known exactly solvable Schrodinger equations are directly related to some 
shape invariant operators, and most of the formulae occurring in the study of these 
quantum systems follow from a small number of mathematical results concerning the 
hypergeometric type operators. It is simpler to study these shape invariant operators 
then the corresponding operators occurring in various applications to quantum mechanics. 
Our systematic study recovers known results in a natural unified way, and allows one to 
extend certain results known in particular cases. 



Nicolae Cotfas / Central European Science Journals () 1-13 



13 



Acknowledgment 

This research was supported by the grant CEx06- 



References 

[1] F. Cooper, A. Khare, U. Sukhatme: " Supersymmetry and quantum mechanics" 
Phys. Rep., Vol. 251, (1995), pp. 267-385. 

[2] N. Cotfas: "Shape invariance, raising and lowering operators in hypergeometric type 
equations", J. Phys. A: Math. Gen., Vol. 35, (2002), pp. 9355-9365. 

[3] N. Cotfas: "Systems of orthogonal polynomials defined by hypergeometric type 
equations with application to quantum mechanics", CEJP, Vol. 2, (2004), pp. 456- 
466. See also http://fpcm5.fizica.unibuc.ro/~ncotfas . 

[4] J.W. Dabrowska, A. Khare, U. Sukhatme: "Explicit wavefunctions for shape- 
invariant potentials by operator techniques" , J. Phys. A: Math. Gen. , Vol. 21, (1988), 
pp. L195-L200. 

[5] L. Infeld and T.E. Hull: "The factorization method ", Rev. Mod. Phys., Vol. 23, 
(1951), pp. 21-68. 

[6] M.A. Jafarizadeh and H. Fakhri: " Parasupersymmetry and shape invariance in 
differential equations of mathematical physics and quantum mechanics", Ann. Phys., 
NY, Vol. 262, (1998), pp. 260-276. 

[7] A.F. Nikiforov, S.K. Suslov, V.B. Uvarov: Classical Orthogonal Polynomials of a 
Discrete Variable, Springer, Berlin, 1991. 



